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Abstract. Using the complete numerical solution of a time-dependent three-dimensional mean-field model 
we study the Josephson oscillation of a superfluid Fermi gas (SFG) at zero temperature formed in a 
combined axially-symmetric harmonic plus one-dimensional periodic optical-lattice (OL) potentials after 
displacing the harmonic trap along the axial OL axis. We study the dependence of Josephson frequency 
on the strength of the OL potential. The Josephson frequency decreases with increasing strength as found 
in the experiment of Cataliotti et al. [Science 293 (2001) 843] for a Bose-Einstein condensate and of the 
experiment of Pezze et al. [Phys. Rev. Lett. 93 (2004) 120401] for an ideal Fermi gas. We demonstrate 
a breakdown of Josephson oscillation in the SFG for a large displacement of the harmonic trap. These 
features of Josephson oscillation of a SFG can be tested experimentally. 

PACS. 03.75.Ss Degenerate Fermi gases, 03.75.Lm Tunneling, Josephson effect, Bose-Einstein conden- 
sates in periodic potentials, solitons, vortices, and topological excitations, 03.75.Kk Dynamic properties of 
condensates; collective and hydrodynamic excitations, superfluid flow 



1 Introduction 



The observation of an oscillating Josephson current across 
an one-dimensional (ID) periodic array of potential wells, 
usually generated by a polarized standing-wave laser field 
and commonly known as an optical-lattice (OL) poten- 
tial, in a trapped "cigar-shaped" Bose-Einstein conden- 
sate (BEC) by Cataliotti et al. Q] was the first manifesta- 
tion of this phenomenon in trapped neutral bosons. Until 
then the Josephson effect was confirmed in superconduc- 
tors with charged electrons and in liquid helium [Ull]- In 
the experiment of Cataliotti et al. [1] a Josephson oscil- 
lation was initiated in a re pulsive 87 Rb BEC formed m a 
ID periodic OL plus an axially-symmetric harmonic po- 
tentials by suddenly displacing the harmonic trap along 
the axial direction. They found that the Josephson fre- 
quency reduced with the increase of the strength of the 
OL potential. Also a breakdown of Josephson oscillation 
was observed for large displacement of the OL potential. 

There have been theoretical studies of Josephson oscil- 
lation using the numerical solution of the time-dependent 
mean-field Gross-Pitaevskii (GP) equation [3] in one (ID) 
[1,4! and three (3D) [5l[6] dimensions to understand differ- 
ent experimental features pQ. There have also been other 
theoretical studies of Josephson oscillation in a trapped 
BEC [Il[8] using different approaches and under different 
conditions distinct from the experiment of Cataliotti et al. 
[I]. Moreover, there have been many interesting studies on 
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different aspects of Josephson oscillation in a cold Fermi 
gas under diverse conditions of trapping [9] . 

Recently, in another experiment on ideal Fermi gas of 
K atoms in an axially-symmetric plus OL trap Pezze 
et al. [10] observed the oscillation of the Fermi gas after 
a sudden displacement of its equilibrium position along 
the lattice for various OL strength. They measured the 
frequency of oscillation for various OL strength. They also 
explained the result of experiment using a semiclassical 
theory. 

Due to a strong repulsive Pauli-blocking interaction at 
low energies among spin-polarized fermions, there cannot 
be an evaporative cooling leading to a quantum degenerate 
Fermi gas (DFG) [11]. Trapped DFG has been achieved 
only by sympathetic cooling in the presence of a second 
boson or fermion component. Recently, there have been 
successful observation [Hl[l^[TMT4] and associated exper- 
imental [mOIllIll and theoretical [20lJ2T][22l[T8|[T9] studies 
of BEC-DFG mixtures. Different experimental groups [TT1 
[T2 i rr3 i ri4] observed the formation of a DFG in the follow- 
ing systems: 6 < 7 Li [13], 23 Na- 6 Li [H] and 87 Rb- 40 K [HI 
ITS] . Also, there have been studies of a degenerate mixture 
of two hyperfine-spin components of fermionic 40 K |llj 
and 6 Li [T^] atoms. Later on, the observation [53] of the 
transition of a DFG to a Bardeen-Coopcr-Schrieffcr (BCS) 
superfluid fermion gas (SFG) of opposite hyperfine-spin 
orientation by manipulating the fermion-fermion interac- 
tion using a Fcshbach resonance has opened the possibil- 
ity of controlled study of a BCS superfluid formed under 
the action of a weak fermion-fermion attraction. The use 
of a Fcshbach resonance has also allowed the observation 
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of BCS-BOSE crossover [24] in two-hyperfinc-componcnt 
fermion vapors of 40 K [2tp3] and 6 Li [26] atoms. 

There have been theoretical models of a DFG using a 
mean-field hydrodynamic approch based essentially on the 
Thomas- Fermi- Weizsacker approximation [27lfT8] . Such mod- 
els have provided results [19] for collapse in a boson-fermion 
mixture of 87 Rb- 40 K in qualitative agreement with the ob- 
servation by Modugno et al. [T5][22] . There have also been 
studies of fermionic bright [55] , dark [35] , and gap solitons 
|30] in a boson-fermion mixture where the bosonic compo- 
nent is treated by the mean-field GP equation [3] and the 
fermionic component by a hydrodynamic model [191118] . 
(The Fcshbach resonances in the boson-fermion systems 
23 Na- 6 Li and 87 Rb- 40 K have been observed experimen- 
tally |31j and can be used in a controlled experiment of 
solitons in such a mixture.) The results of such studies 
are in agreement with ab initio studies of gap [32j and 
bright [33| solitons based on properly antisymmetrized 
many-body approach. 

Although a mcan-ficld-hydrodynamical approach to a 
DFG has its limitations (as in the absence of a coher- 
ent phase it may determine only the fermion density), 
a Ginzburg-Landau-type mean-field description [34j of a 
BCS superfluid based on a Lagrangian density is theoret- 
ically well founded and widely appreciated and has been 
used [301135] in dealing with superfluid boson-fermion mix- 
tures (such a description for a SFG determines the fermion 
density as well as the phase of a complex coherent or- 
der parameter). A rigorous antisymmetrized many-body 
approach for a SFG becomes unfeasible as the number 
of fermions increases, where the simplified mean-field ap- 
proach of a BCS superfluid in terms of a one-body wave 
function of Cooper-paired fermions is of advantage. 

In this paper we study the Josephson oscillation of a 
BCS superfluid at zero temperature in an OL potential 
using the complete numerical solution of a 3D mean-field 
model. As in the experiment on Josephson oscillation of 
a BEC [15], we consider a SFG formed in a cigar-shaped 
axially-symmetric parabolic trap with an added OL poten- 
tial along the axial direction. The Josephson oscillation is 
initiated by giving a translation of the parabolic trap in 
the axial direction upon the formation of the SFG. The 
SFG acquires energy in the process to initiate the Joseph- 
son oscillation. We study the variation of the frequency of 
Josephson oscillation for different OL strength. We also 
study the breakdown of smooth oscillation for large ini- 
tial displacement of the parabolic trap. We compare our 
results for frequency with experimental results on super- 
fluid Bose 87 Rb atoms p] as well as on Fermi 40 K atoms 



The present model is derived as the Euler-Lagrange 
equation of a Lagrangian density using the well-known 
energy density of a SFG [36]. The Euler-Lagrange equa- 
tion so obtained is a nonlinear partial differential equation 
with a nonlinearity of power 7/3. Such an equation has 
also been used for a DFG with a different coefficient mul- 
tiplying the nonlinearity (the coefficient in the two cases is 
different due to the Cooper pairing in a BCS superfluid). 



In Sec. II we present the mean-field model for a SFG 
mixture. In Sec. Ill we present numerical results for Joseph- 
son oscillation of a SFG in a combined harmonic plus ID 
OL potentials as a function of the strength of the OL when 
the harmonic trap is displaced through a small distance 
along the optical axis. We also illustrate the breakdown 
of Josephson oscillation in a SFG when the oscillation is 
initiated by giving a large displacement of the harmonic 
trap along the optical axis. Finally, in Sec. IV we present 
a summary of our study. 



2 Mean-field Model for Fermi superfluid 

Let us consider a BCS superfluid Fermi gas (SFG) of N 
Cooper-paired fermions of mass m. The energy density 
of the system is given by £ = 3aep/5 [36], where ep = 
h 2 k F /(2m) is the Fermi energy with kp the Fermi mo- 
mentum and a the atomic density. (Here we neglect small 
corrections to this expression due to the residual fcrmion- 
fcrmion interaction usually expressed as an expansion in 
scattering length.) Modification to this expression appro- 
priate to study a BCS-Bosc crossover [24] has also been 
suggested [37] • The atomic density for the superfluid corre- 
sponding to the isotropic distribution (spherical Fermi sur- 

face) is a = 2(2^)~ 3 / fcF 4nk 2 dk = (Sir 2 )" 1 (2me F /h 2 ) 3/ \ 
where the effect of Cooper pairing is included in the ex- 
tra factor of 2. This relation can be inverted to yield 



ep = h (37r a) ' 3 /(2m), which leads to the following ex- 
pression for energy density of the system 
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In the spirit of the Ginzburg-Landau theory [51] the 
SFG can be described by a complex order parameter 
such that a = 'P 2 . In terms of this order parameter, the 
Lagrangian density of the superfluid can be written as 
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where V(r) is the external potential, and m c g is the ef- 
fective mass of superfluid flow as in the Ginzburg-Landau 
theory (The exact value of m e g is not known but there 
is evidence [34] that for a Cooper-paired SFG it is twice 
the fermion mass and we shall use this value in the fol- 
lowing: TO ff = 2m). The Euler-Lagrange equation of La- 
grangian density ([2]) becomes the following 3D nonlinear 
Schrodinger equation with a repulsive nonlinear term of 
power 7/3: 



iW t = 



2m off 



-V 2 f 



2m 



(3tt 



2 ^ 2/3 i.r/i 4 /3 



+ V(r)9 , (3) 



with normalization J |^(r, t)\ 2 dr = N. In the presence 
of the combined axially-symmetric and periodic OL po- 
tentials [l][5] V(r) — \m e su> 2 {p 2 + v 2 z 2 ) + Vq^, where 
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to is the angular frequency of the harmonic potential in 
the radial direction p, voj that in the axial direction z, 
with v the aspect ratio, and V ^i = sEncos 2 (kLz) is the 
OL potential created with the standing-wave laser field of 
wavelength A, with Er = % 2 k\/{2m), = 2n/\, and s 
(< 12) the strength. Although, for small s one should have 
a Joshepson oscillation, for large s values Joshepson oscil- 
lation should terminate, as in the case of a BEC [38], due 
to a superfluid to Mott-insulator transition. (This transi- 
tion cannot be studied with the mean-field equations and 
field-theoretic analysis is needed for its understanding.) 

In the axially-symmetric configuration, in the zero an- 
gular momentum state the fermion order parameter can 
be written as ^(r,i) = ip(p,z,t), where < p < oo and 
— oo < z < oo. Now transforming to dimcnsionlcss vari- 
ables p = \/2p/l, z = \/2z/l, t = tuj, I = \/Ti/ (m c guj), 



and ip(p, z; r) = 

mm 



P/(NV8)i/j(p,z;t), Eq. © becomes 
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where the nonlinearity parameter n = m c g (37t 2 7V) 2 / 3 /to. 
In terms of the ID probability P{z 1 1) = 2ir dp\ip(p, z, t)\ 
the normalization of the wave function is given by dzP( 
1. 

In many of the experiments with DFG [mflMTO] . 40 K 
fermion atoms were used and in this study we also con- 
sider 40 K atoms and take m to be the mass of K atoms. 
As in the experiment of Cataliotti et al. [I] the axial and 
radial trap frequencies are taken as vuj = 2tt x 9 Hz and 
uj = 2tt x 92 Hz, respectively, with v = 9/92 ps 0.1. For 
40 K, the harmonic-oscillator length / = \Jfij (m c ffo;) ps 1 
pm and the present dimcnsionlcss length unit corresponds 
to 1/V2 ps 0.7 /im. The present dimensionless time unit 
corresponds to uj^ 1 = l/(2ir x 92) s = 1.73 ms. Although 
we perform the calculation in dimcnsionlcss units using 
Eq. 0, we present the results in actual physical units 
using these conversion factors consistent with 40 K atoms. 
We take the wavelength A of the standing-wave laser beam 
to make the OL potential to be A = 700 nm, so that 
the dimensionless laser wave length A = V^X/l ~ 1 
and the dimcnsionlcss standing-wave energy parameter 
Ef>/(huj) = 47t 2 /Aq. Hence in dimensionless units the OL 
potential of Eq. J4| is 
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Though most of our calculation was done with above set 
of parameters, for a comparison with the experiment on 
Fermi 40 K atoms of Pezze et al. [TU], we also repeated our 
calculation with trap frequencies vlo = 2ir x 24 Hz and 
uj = 2ir x 275 Hz, and OL wavelength A = 863 nm. 

We solve the three-dimensional Eq. numerically us- 
ing a split-step time-iteration method with the Crank- 
Nicholson discretization scheme described recently [55] . 



We discretize the GP equation typically with time step 
0.001 and space step 0.1 spanning p from to 7 /im and 
z from —120 pm to 120 pjn, although, sometimes we used 
smaller steps for obtaining convergence. Equation (j4]) is 
then solved by time iteration starting with the known har- 
monic oscillator solution for n = 0: ip(p, z) = [v/fin 3 )] 1 / 4 
p e~{P +uz )/ 4 with chemical potential p, = (1 + v/2) 
[40] . For a typical cigar-shaped geometry with v ~ 0.1 
PP, p{~ 1) is much smaller than the typical depth of 
the OL potential wells E r /{Tiuj) = 47r 2 /A§ ~ 40 so that 
p << Er/ (Tiuj) and the passage of fermion atoms from one 
well to other can only proceed through quantum tunneling 
5,6. During the time iteration of Eq. ^ the nonlinearity 
n as well as the OL potential parameter s are slowly in- 
creased by equal amounts in 100000 steps of time iteration 
until the desired value of nonlinearity and OL potentials 
are attained. Then, without changing any parameter, the 
solution so obtained is iterated 50000 times so that a sta- 
ble solution is obtained independent of the initial input 
and time and space steps. The solution then corresponds 
to the bound SFG under the joint action of the harmonic 
and OL potentials. 



3 Numerical Results 

2 F^jrst we consider a SFG formed in the combined harmonic 
andjicriodic OL potentials for a specific nonlinearity n. 
We study the formation of a SFG in the combined har- 
monic and OL potentials of Eq. ([5]) for a range of values 
of s. In Fig. 1 (a) the plot of \ip(p, z)\ vs. p and z illus- 
trates the profile of the SFG for N = 40, n = 225 and 
s = 4. From Fig. 1 (a) we see that the SFG has the shape 
of a cigar (of length 80 pm and transverse radius 3 pm) 
cut into narrow slices with the OL barriers separating the 
slices. The large number of maxima and minima due to 
the OL potential is not clearly visible in this plot. The 
maxima and minima in the axial direction are clear in the 
plot of \ip(p, z)\ vs. p and z in Fig. 1 (b) where we show the 
central part of function \tp(p, z)\ for 3 pm > z > —3 pm. 
In this interval of z, there are about 16 wells of the OL 
potential and as many maxima and minima in \ip{p, z)\. 

Now we consider an oscillating SFG in the combined 
harmonic and periodic OL potentials. If we suddenly dis- 
place the harmonic trap along the lattice axis by a small 
distance after the formation of the SFG in the combined 
potentials, the SFG will acquire a potential energy, be out 
of equilibrium and start to oscillate. As the height of the 
potential-well barriers of the OL potential is much larger 
than the chemical potential of the system, the atoms in the 
SFG will move by tunneling through the potential barri- 
ers. This fluctuating transfer of atoms across the potential 
barriers is due to Josephson effect in a neutral quantum 
SFG. 

With the SFG of Fig. 1 (a) we next study its Joseph- 
son oscillation when the harmonic potential is suddenly 
displaced along the axial direction by 14 pm. The SFG 
now acquires an added potential energy which it can use 
to execute a Josephson oscillation along the axial direc- 
tion. The Josephson oscillation is best studied numerically 
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Fig. 1. (a) The profile of the SFG function \ip(p, z)\ vs. p and 
z formed in the joint optical-lattice and harmonic potentials 
for N = 40, n = 225 and s = 4. (b) The profile of the SFG 
function \ip(p, z)\ of (a) vs. p and z for 3 pm > z > —3 pm, 
which clearly exhibits the maxima and minima of the SFG 
function along the axial direction. 



from the contour plot of the ID probability P(z,t) vs z 
and t exhibited in Figs. 2. These plots clearly show the 
central position of the SFG along the axial z direction. 
In the present simulation we take different values of the 
OL strength s. These contour plots are very useful to find 
the Josephson frequencies. From Figs. 2 the periods of 
Josephson oscillation are easily read off and the frequen- 
cies calculated for different OL strengths s. 

In Fig. 3 (a) we plot the Josephson frequencies vs. OL 
strength s. Specifically, in addition to the present calcula- 
tion for SFG we also show (i) the experimental frequencies 
of Cataliotti et al. p] for a repulsive BEC of 87 Rb atoms, 
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Fig. 2. (Color online) The contour plot of the ID probability 
density P(z, t) of the SFG with n = 225 executing a Josephson 
oscillation when the harmonic trap is suddenly displaced along 
the axial direction through a distance 14 pm at t = for 
OL strength s = (a) 6 and (b) 8. The period of Josephson 
oscillation can be obtained from these plots. 



(ii) the 3D simulation for a BEC from Ref. [6], and (ii) 
the hydrodynamical calculation for a BEC from Ref. [TJ. 
In Fig. 3 (b) we show the Josephson frequencies vs. OL 
strength with the parameters of the experiment of Pezze et 
al. [10] for an ideal Fermi gas of 40 K atoms. In this figure 
we compare our results with experiment and a semiclassi- 
cal theoretical calculation [TU] . 

The most interesting conclusion from Figs. 3 is that 
the present frequencies of 3D simulation for Josephson 
oscillation of a SFG are practically the same as the fre- 
quencies for Josephson oscillation of a BEC and close to 
those for oscillation of an ideal Fermi gas. For s — the 
OL potential is absent and the BEC and the SFG execute 
free oscillation with the frequency of the axial potential. 
From Figs. 3 we find that the Josephson frequency reduces 
with increasing OL strength. As the Josephson oscillation 
takes place by quantum tunneling of the atoms through 
the OL barriers, this oscillation is bound to be reduced as 
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Fig. 3. (Color online) The frequency of the atomic current in 
the array of Josephson junctions as a function of OL strength 
s compared with (a) experiment on boson S7 Rb and other the- 
ories (b) on fermion 40 K and other theory. In (a) solid square 

— present result for SFG; solid circle — experiment of Catal- 
iotti et al. [I] for BEC with repulsive 87 Rb atoms; full line 

— mean-field simulation for BEC [B]; dotted line — hydrody- 
namical calculation for BEC [7]. In (b) solid square — present 
result for SFG; solid circle — experiment of Pezze et al. [ID] 
for fermions with 40 K atoms; full line — semiclassical theory of 
Pezze et al. [TD] for fermions; dotted line — joins present points 
to guide eye. 



the height of the OL barriers is increased. The reduction 
of Josephson oscillation with the increase of the parame- 
ter s results in a reduction of Josephson frequency in Fig. 
3. The agreement of the frequencies of present model for 
superfluid Fermi gas with the experiment on ideal (non- 
supcrfluid) Fermi gas is remarkable. A future experiment 
on a superfluid Fermi gas may reveal subtle differences 
in behavior, if any, of such oscillations in superfluid and 
nonsuperfluid Fermi gas. 
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Fig. 4. (Color online) The contour plot of the ID probability 
density P(z, t) of the SFG with n — 225 released after hold 
times Th = (a) ms, (b) 10 ms, (c) 30 ms, and (d) 50 ms 
when the harmonic trap is suddenly displaced along the axial 
direction through a distance 100 /im for OL strength s = 8. The 
three clean interference trails in (a) demonstrates the phase 
correlation in the initial SFG. With the increase of hold times 
in the displaced trap the clean interference trail is slowly lost 
in (b) and (c) signalling a loss in phase correlation. The phase 
correlation in the SFG is completely lost in (d) after a hold 
time of 50 ms. 



The SFG formed on the combined OL and harmonic 
traps considered so far in Figs. 1 and 2 is a phase-correlated 
quantum fluid with the atoms freely moving by quantum 
tunneling from one OL site to another at zero tempera- 
ture and one has a phase-coherent macroscopic state. This 
will lead to an interference pattern when the SFG is re- 
leased from the OL trap. However, the details of the pat- 
tern may change from one experiment to another. When 
the harmonic trap is displaced through a small distance 
as in Figs. 2, the SFG executes Josephson oscillation by 
quantum tunneling. 

Next we consider an account of the breakdown of Joseph- 
son oscillation in a SFG when the harmonic potential is 
displaced by a very large distance (~ 100 /im) along the 
periodic OL potential. In that case for small times (t < 300 
ms) no Josephson oscillation of the type observed in Figs. 
2 was found. The SFG was found to remain virtually fixed 
at z = for small times. At very large times it moved 
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very slowly to the new center of equilibrium. However, the 
phase correlation of the SFG is lost after staying a signif- 
icant time in the displaced trap. To demonstrate the de- 
struction of phase correlation in the displaced trap we rely 
on the disappearance of the interference pattern formed 
upon the release of the SFG from the combined harmonic 
and OL traps after a hold time of 1\ in the displaced trap. 

The formation of a BEC on a combined harmonic plus 
OL potentials and the interference pattern formed upon 
the release of this BEC from the confining traps have been 
studied both theoretically [4][42j|4l] and experimentally p] 
38,43]- We present a similar study for the SFG. The pe- 
riodic OL potential cuts the SFG into several pieces at 
different sites maintaining phase correlation among them. 
As a result when the SFG is released from the traps it 
expands freely and a matter-wave interference pattern is 
formed in a few milliseconds. The atom cloud released 
from one lattice site expands, overlaps and interferes with 
atom clouds from neighboring sites to form a robust in- 
terference pattern, consisting of a central peak and two 
smaller symmetrically spaced peaks moving in opposite 
directions along the OL axis [Pjl3"8]l43| . Since the lattice 
transfers momentum to the SFG in units of 2pn = 2/i/A, 
the recoil velocity of each the two side peaks is given by 
Vr = 2 PR /{2m) = 2h/(2mX) PEE]. Using I 2 = %/{2ww), 
we have vr = inPoj/X w 10 mm/s, where we used the 
numerical values I ~ 1 /im, u> = 1.73 ms, and A « l/y/2 
/.im. 

To demonstrate the loss of phase correlation after dis- 
placing the harmonic trap along the OL direction through 
100 ^m, in Figs. 4 we exhibit the contour plot of the ID 
probability density P(z,t) of the SFG with n = 225 re- 
leased after different hold times T% in the harmonic plus 
OL traps. In Fig. 4 (a) we exhibit the result of simulation 
for Th = 0, which has no effect on the phase correlation. 
Upon release from the traps a robust interference pattern 
of the central and two side peaks can be seen in this figure. 
The side peaks move with velocities 10 mm/s and in 12 
ms they move about 120 ^m each as can be seen from this 
figure. The successive figures 4 (b), (c), and (d) show the 
situation for hold times 10 ms, 30 ms, and 50 ms, respec- 
tively. The clean interference pattern of a central and two 
symmetrically moving side peaks is destroyed slowly with 
the increase of hold time. For an intermediate Th (= 10 ms 
and 30 ms) there is partial destruction of the interference 
pattern. For Th = 50 ms, there is no interference pattern 
and the expanding SFG occupies the full region between 
the two side peaks. This simulation illustrates that when 
the harmonic trap is displaced by a large distance, there 
is no Josephson oscillation. 

The breakdown of Josephson oscillation above is due 
to a classical transition from a phase-correlated superfluid 
to an "insulator" resulting in a modulational instability as 
in a BEC [44] . Other mechanisms for the loss of phase co- 
herence in a BEC have also been studied [IS]. The loss 
of phase coherence considered in all these investigations 
[44"]l4"5] originated from a classical superfluid to insulator 
transition, different from a quantum transition of a super- 
fluid to a Mott insulator observed in Rcf . [38] . The classi- 



cal phase transition involves energy, whereas the quantum 
phase transition does not involve a supply of energy. 

A recent investigation [15] indicates that, for a BEC, 
the disruption of the large amplitude oscillations is caused 
by the onset of dynamical instability occurring when the 
quasi-momentum surpasses a critical value towards the 
edge of the first Brilluoin zone. In this regard, it would 
have been interesting to study what is the critical quasi- 
momentum as a function of the nonlinearity of a SFG. In 
another experimental study on superfluid Fermi gas [47] 
with an OL moving at a constant speed, the breakdown 
of superfluidity was measured as function of speed over 
the entire Bose to BCS crossover. A careful theoretical 
investigation on these issues, although beyond the scope 
of present work, are welcome in the future. 



4 Summary and Conclusion 



We performed 3D numerical simulation based on a time- 
dependent mean-field equation to study Josephson oscilla- 
tion of a BCS superfluid in a combined axially-symmetric 
harmonic and periodic OL potentials. The OL potential 
is aligned along the axial direction and is created by a 
standing-wave laser beam. The Josephson oscillation is ini- 
tiated by displacing the harmonic potential along the OL 
axis by a small distance (~ 15 ^m). We study the varia- 
tion of Josephson frequency with the strength of the OL 
potential and find that the frequency decreases with the 
increasing OL strength s. As s increases from 1 to 10 the 
frequency of our SFG model is in agreement with an ex- 
pcricmnt on BEC and on ideal Fermi gas. However, for a 
large displacement (~ 100 fim) of the harmonic potential 
along the OL axis we demonstrate a breakdown of Joseph- 
son oscillation. The SFG then returns very slowly to the 
new mean position without executing Josephson oscilla- 
tion. The breakdown of Josephson oscillation also leads to 
a destruction of phase correlation and consequent loss of 
superfluidity in the Fermi gas and is demonstrated by al- 
lowing the SFG to undergo a free expansion. The absence 
of a distinct interference pattern upon free expansion sig- 
nals the loss of phase correlation due to a classical phase 
transition to an insulating state as opposed to a quan- 
tum superfluid to a Mott insulator transition observed in 
a BEC [38] . Similar superfluid to insulator classical phase 
transition has been observed [Tll4"l] and studied [5ll4"5] in a 
BEC. These features of Josephson oscillation could be ver- 
ified experimentally and will provide a test for the mean- 
field model for a SFG. 

A proper treatment of SFG should be performed us- 
ing a fully antisymmctrized many-body Slater determi- 
nant wave function [20]l33] as in the case of atomic and 
molecular scattering [35]. However, in view of the success 
of the hydrodynamic model in other contexts [I9ll28ll29] 
we do not believe that the present study on Josephson os- 
cillation of a DFG in an OL potential to be so peculiar as 
to have no general validity. 
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